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CNj ' The renormalization method which is a type of perturbation method is extended to a 

tool to study weakly nonlinear time-delay systems. For systems with order-one delay, we 
show that the renormalization method leads to reduced systems without delay. For systems 
with order-one and large-delay, we propose an extended renormalization method which leads 
to reduced systems with delay. In some examples, the validities of our perturbative results 
' are confirmed analytically and numerically. We also compare our reduced equations with 

reduced ones obtained by another perturbation method. 
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§1. Introduction 



Many nonlinear dynamical systems in various scientific disciplines are influenced 
by the finite propagation time of signals in feedback loops. A typical physical system 
is provided by a laser system where the output light is reflected and fed back to the 
cavity. Time delays also occur in other situations. For example, in a traffic flow 
I model including a driver's reaction time,^^ in biology due to physiological control 

■ mechanisms, or in economy where the finite velocity of information processing has 
, to be taken into account. Furthermore, realistic models in population dynamics 

CN I or in ecology include the duration for the replacement of the resources. In some 

■ situations, such as lasers and electro-mechanical systems,^) systems with large-delay 
, appear. For this reason, we need to develop a mathematical tool to study them, 

especially for weakly nonlinear systems as a first step. The main difficulty peculiar 
to systems with delay is its dimensionality. Due to a delayed arraignment in a given 
system, x{t — r) = exp(—rd/dt)x(t), the dimension of the phase space is high. 
^ ' Suppose we add a perturbation term to a given system, the system is not guar- 

■ anteed to be structurally stable. So the perturbation result is, if computed naively, 
plagued with singularities such as secular terms. It has been recognized that these 
singularities in the result of the naive perturbation method can be renormalized away 
by the modification (renormalization) of parameters associated with the unperturbed 
system.^) The modified parameters are governed by the renormalization equations 
that turn out to be slow-motion equations (reduced equations). It is important that 
the prescription of the method does not depend on the details of the system under 
study. To obtain a more useful and sophisticated tool to study weakly nonlinear 
systems, reformulated versions^''^"^^^ of the original method^^ have been proposed. 
It is noted that there are a variety of applications of renormalization methods to 
physical systems, such as plasma physics,^'''' general relativity^*^) and quantum op- 
tics, ^^-^ in addition to studies in standard nonlinear dynamical systems. Although 
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the reformulated version of the renormahzation method that we employ here is easily 
applied to non-chaotic systems^^^'^*^^ and chaotic maps,^^^ we do not know whether 
or not the renormahzation method can be applied to systems with delay. 

The purpose of this paper is to show that the reformulated renormahzation 
method can be applied to weakly nonlinear systems with delay. For systems with 
order-one delay, the method leads to reduced systems without delay. For systems 
with large-delay, the method should be extended and the application of the extended 
method leads to reduced systems with delay. Our extended method can also be 
applied to systems with order-one delay, the resultant reduced equations are different 
from those obtained by the use of the conventional renormahzation method. As 
mentioned, a time-delay term of a given system makes the dimension of the phase 
space high. Even in such a case, our method can lead to reduced equations. 

The organization of this paper is as follows. In the next section(^, We show 
that our conventional renormahzation method can lead to reduced equations for 
systems with order-one delay. The validities of our analyzes are shown. In ^ we 
propose an extended version of our reformulated renormahzation method so that we 
deal with a system with large-delay. The definition of large-delay is to be given in the 
beginning of the section. We show that the extended method can also be applied to 
systems with order-one delay. The validity of the extended method is also discussed. 
Finally, in ^ we discuss the features of our methods and conclude our study. 



§2. Conventional Renormalization Method 



In this section, using the conventional renormalization method, we first analyze 
a linear system that has an oscillatory solution, and show that our perturbative 
analysis is in agreement with the exact solution analytically. By the conventional 
method we mean the method proposed in Ref.16). Next, we study some classes of 
weakly nonlinear systems using our renormalization method. Our classes include 
a nonlinear oscillator, a laser model, systems with many degrees of freedom, and 
spatially extended systems. In some examples, we show that our analyzes are valid 
by comparing with the numerical simulation or the previous studies. 

2.1. Linear system 

To show that our renormalization method^^^ can be applied to systems with 
delay, we consider the following system as an example, 

^^^4P-+uj^x{t)+ex{t-r) =0, (2-1) 

where uj{g M) is a parameter, r(G M) represents the time-delay, and e(G M) is a 
small parameter {\e\ <C 1). In this system, there is an oscillatory solution that is 
analytically expressible without any approximation. The exact solution is written as 



x{t) = Aexp(it\/u;2 - e) + c.c, (2-2) 

under the condition 



" :. (2-3) 
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Here c.c. represents for the complex conjugate terms of the preceding expression, 
A(g C) is the integration constant. 

Let us derive a perturbation solution of (j2-ip using our renormalization method. 
In this perturbative analysis, we do not use the exact solution (j2-2p . As well as in 
the case of a differential equation without delay, we first find the naive perturbation 
solution, x{t) = x^'^\t)+ex'^^\t)+e'^x^'^\t) + 0{e^). This naive perturbation solution 
is obtained by solving the following equations, 

Lx(°)(t) = 0, Lx^^\t) = -x(^-i)(t-r),(j = 1,2,...) 
Lxit) := (^^+^2^x(t). 

The solutions are obtained as following 

= Ae^'^* + c.c., 
^(i)(t) = !^ite^-(*-)+c.c., 

^(2) (i) = Z4 f t2 _ 2rt + U) e^-^*-^'-) + c.c. , 

where A{£ C) is the integration constant of the solution of the unperturbed system, 
x^^\t). Note that the solutions x^^\t),{j > 1) contain the terms const. exp(ia;t) 
and const.exp(— icjt). We assume that these terms are included in Aexp{iLot) and its 
complex conjugate term in x^'^^t). Apparently, the validity of the naive perturbation 
solution is invalid in the regime t > 0{l/e), due to the secular terms (oc et, cx e^t^ 
etc.). 

The renormalization method removes the secular behavior in a systematic way. 
We define the renormalized variable A{t) up to 0{e'^), 

A{t) ■.= A + e—te-'^'' + e^Z^ (t"^ - 2rt + -t] e'^^'^^ (2-4) 
2lo ou;^ \ io J 

Note that this definition is a form of a near-identity transformation at the constant 
A,^^^ and that the naive perturbation solution is expressed in terms of the renormal- 
ized variable, ^ 

x{t) = A{t) ex.p{iujt) + c.c. + 0{e^). (2-5) 

We construct the equation which A{t) should satisfy perturbatively. Such an equa- 
tion is our renormalization equation. From Eq. (j2-4p . we obtain the following two 
relations, 

A{t + a) - A{t) = e—ae-'^' 
2uj 

+e^^(2ta + a'^-2ra + -a\-'^'^' + 0{e^). (2-6) 



and 



yl = ^(t) -e^^ie-^'^" + 0(e2), (2-7) 
2uo 



4 



S. Goto 



where cr(G M) is a parameter. Substituting Eq. (l2-7p into Eq. (l2-6p . we have the 
approximate closed relation 

The renormalization equation is obtained in the limit o" ^ as 

+'^^*\4^"&7j' • ^^-^^ 

The solution of Eq. (j2-8p is given by 

= I(0)e'^W, (2-9) 

To compare the solution of our renormalization method with the exact solution 
(j2-2p . we restrict ourselves to the approximate solution imposed on the condition 
(IMl) . Using Eq. ([231), we can rewrite Eq. ([23]) as 

Ait) = 1(0) exp + + • (2-10) 

In terms of x{t), we obtain the approximate solution using Eqs. (j2-5p and (j2T0p . 

= 1(0) exp jii - ^ - ^ + 0(e^)) | + c.c. + 0{e^). (2-11) 
In fact, due to the relation 

Eq. (|2TT]1 is the same as Eq. ([22]) up to ©(e^). 

2.2. Nonlinear single oscillator 

Let us consider a nonlinear equation including a time-delay term, which is to 
show that our renormalization method can be applied to such a system. The system 
which we study here is 



dt 



+ ax{t) + f3x{t -r)= e{-iix{t) - -fsx^t)), (2-12) 



where a,/9,7i(G M) and 73 (G M) are parameters. The value of r(E M) represents 
the time-delay, and e(G M) is the small parameter. It is noted here that the unper- 
turbed system has an analytically expressible oscillatory solution when the following 
condition is satisfied, 

^="7,f°f ■ (0'>^% (243) 
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We restrict ourselves to the case that this condition is satisfied. The oscillatory 
solution of the unperturbed system is given by 



X 



=^e^'^* + c.c., 



where 00 := \l fP' — a^, A(s C) is the integration constant, and the relation iw + a = 
— /3exp(— iwr) is satisfied. Although we can analytically express the exact solution 
of the unperturbed system, it is difficult to express the exact solution analytically in 
the case e 7^ 0. To investigate the effect of the perturbation term, we construct the 
renormalization equation for Eq. (j2-12p . The procedure to obtain the renormalization 
equation is the same as that described in ^2.11 

The naive perturbation solution, x{i) = x^^'^it) + ex(^)(t) + ©(e^), is obtained 
by solving the following equations, 

L,x(°)(t) = 0, Lrx^^\t) = 7ix(°)(t) - 73x(°)3(t), 
Lrx{t) := ^^^^ + ax{t) + (ix{t - r). 
The solutions are given by 

=^e*^* + c.c., 
x«(t)=t21^i^,^I^e- + c.c.. 

The naive perturbation solution includes the secular term. To remove the secular 
term, we define the renormalized variable j4(t). 

^ ' l + r(a + iu;) 

The renormalization equation is obtained as 

dA{t) _^ 7iI(t)-373|I(t)pI(t) 
dt 1 + r(a + iio) 

This system can be split into two parts: dynamics described by its amplitude and 
phase, 

-?>i,QrAr^-^\. (2-14) 



dt V 373 

d^ ^ 3£73u; / 2 _ 71 \ 



,0 . - — ■ (2T5) 

dt (1 +ra)2 + (ra;)2 V 873; ^ ' 

Here, A{t) := ii(t)e'<^(*) and Q := 2e(l + ra)/{ {l + ra)^ + {rujf}. It turns out that, 
from Eqs. (|2-14p and (j2-15p . R = R^, (i?* := 71/(873)) is a stable fixed point when 
71/73 > and 73Q > 0. It is thus expected that the limit-cycle oscillation appears 
when these conditions are satisfied. The amplitude of this limit-cycle oscillation is 
given by 2R^: due to the relation, x{t) ^ 2R{t) cos(u;t -|- (l){t)). Fig{T] shows that our 
analysis is valid. 
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Fig. 1. The time sequence of the delay system (|2-12|) . The values of the parameters are a = 
0.1, P = 0.2, 7i = 3, 73 = 1 and e = 0.01. The value of the time-delay r is set to that of (|2-13p . 
The initial condition is x(t) — 0.2 for — r < t < 0. The amplitude of the limit-cycle is 2, as 
predicted for this case [ 2_R« = 21^71/(873) — 2, see text ]. The numerical simulation was done 
using a fixed-step fourth-order Runge-Kutta method with linear interpolation for the required 
two midpoint evaluations of the delayed variable. 



2.3. Lang-Kobayashi phase equation 

In this subsection, we show that our analysis can lead to a set of reduced equa- 
tions which describe dynamics of semiconductor lasers with feedback. In Ref.22), 
they have analyzed the following delay equation 

+ Lo^^" +<P' -A + ioAi cos[^(5 - 0) - ^{S)] = 0, (2-16) 

where ^,A,Ai{g M) are parameters, uj{£ M) is the small parameter, and the prime 
means differentiation with respect to S. This equation is obtained from the Lang- 
Kobayashi equations in the following conditions:^^) the small ratio of the phonon and 
carrier lifetimes and the relatively large value of the linewidth enhancement factor. 
We derive a reduced equation from Eq. (|2-16p using our method and compare our 
result with that obtained by using the multiple-scale method in Ref.22). 

First, the naive perturbation solution, ^{S) = ^^^\S) w^(^)(S') + ©(w^), is 
obtained by solving the following equations 

L^{0)(5) = A, (2-17) 
= -^^^^\S)-Aicos[^'^°\S-e)-<P^^\S)], (2-18) 

The solution to the unperturbed system (j2-17p . 'P^^\S), is obtained as 

<p(0){S) = :^e*(^+^) + c.c. + SA + B, (2-19) 

where A,v{g M), and M) are the integration constants. Substituting Eq. (j2T9p 
into Eq. (j2-18p . we obtain the following equation 



= -CAcos{S + v)- AiJo{D) cos{OA) 
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J.n{-e/2+S+v) 



-yli|cos(0Z\) Jn{Dy 

in{-0/2+S+v) ^ 

+ sm(0Z\) ^ Jn{D) : + c.c. L (2-20) 

n=l,3,-- ^ ^ 

Here D is defined by D := 2^sin(6'/2), J„ denote the n-th order Bessel functions, 
and we have used the fohowing relation in deriving Eq. (j2-20p . 

e'^^^^ = jQ{z) + 2i J„(z)sin(n0) + 2 ^ J„(z) cos(n0). 

n=l,3,--- n=2,4,-- 

The solution of Eq. (j2-18p is given by 

= -^^cos{S + v)- AiMD)Scos{eA) 
+SAiJi{D) sin(6)Z\) sin(-0/2 + S + v) 

Jn{-0/2+S+v) 



^|cos(0Z^) ■ n 2^ 

I ^-^ mil — 

n=2,4,--- ^ ^ 

m(-e/2+5'+-u) ^ 

-sMOA) Y Jn(D) +C.C. . (2-21) 



n=3,5,- 

This naive perturbation solution includes the secular terms (oc ujS). 

Next, we remove these secular terms using our renormalization method. We 
define the renormalized variables C{S){£ C) and B{S){£ M) as follows 

C{S) := A-'^Ua + 2iAiJi{D{A)) sin(0Z\)e-^®/2 



2 

B{S) -.= 3- ujSAiJo{D{A)) cos(0Z\). 



The set of the renormalization equations up to 0{uj) is obtained as 

= -1^iC{S) - iu:AMD{\C{S)\)) sin(0Z^)e-^®/2, (2-22) 
-ujAiJQ{D{\C{S)\))cos{eA). (2-23) 



dS 

Here, we compare the renormalization equations (|2-22p and (I2-23P with the 
reduced ones obtained by the multiple-scale method. Using the decomposition 
C{S) = A{Sy^^^\ {A{S) e ^,v{S) G R), we obtain 

= -^^1(5) - LoA^ sm{eA)MD{A{S))) sin(0/2), (2-24) 

^ = sm{0A)MD{A{S))) cos(0/2). (2-25) 

dS A{S) 
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When we introduce the slow variable C := ujS, the renormalization equations (l2-23p . (|2-24p 
and (j2-25p become the reduced equations derived in Ref.22). This comparison shows 
that our analysis is consistent with one by a traditional perturbation method, and 
that our method can lead to the reduced equations from a physical system. 

2.4. Weakly nonlinear lattice 

In this subsection, we show that our method leads to a discrete complex Ginzburg- 
Landau equation from a weakly nonlinear lattice with delay. In this lattice system, 
the finite propagation time of motion from the nearest oscillators is taken into ac- 
count. Studying the derived reduced system, we predict the stability of a trivial 
solution, and this prediction is confirmed numerically. 

The weakly nonlinear oscillator which we study here is given by 



dxjjt) 

dt 
dpjjt) 

dt 



: -f2'^Xj{t) 

+£i^u(^Xj+i{t 



r) + Xj-i{t — r) — 2xj{t) 



ax^j{t) 



(2-26) 



(2-27) 



where a, M) are parameters, e(G M) is the small parameter, and r(G M) represents 
the time-delay. The variables Xj{t){G M) and pj{t){G M.) denote the displacement 
and momentum of the single oscillator located at lattice site ^) respectively. It 
is noted that this given system becomes a Hamiltonian system when r = 0. Using 
the conventional renormalization method, we derive the reduced system here. 

First, the naive perturbation solutions Xj{t) = x^^\t) + £x^^\t) + 0{e'^) are 
obtained as 



+£■ 



te' 



2iQ 



u{e-'"'' {Aj+i + Aj^i) - 2Aj} - 3a\Aj\'^Aj 



+ C.C. (2-28) 



Here Aj{£ C) are the integration constants, and the higher harmonic terms in Xj (t) 
are omitted. 

Second, from Eq. (j2-28p . the renormalized variables are defined as 



Aj{t) := Aj + e 



2iQ 



u{e-^^''{Aj+i + Aj^i) - 2Aj} - 3a\Aj\'^Aj 



From the definitions, we have the relation Xj{t) « Aj{t)e'^^ + c.c, and the renor- 
malization equations 

dAj{t) _ . 



dt 



2if2 



Ke-^^'^(^,+i(t) + - 2A,{t)} 



-3a\Aj{t)\'^Aj{t) 



(2-29) 
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The system (12-29I) becomes the discrete nonhnear Shcrodinger equation, a Hamil- 
tonian system, in the case r = 0. When r 7^ 0, Eq. (j2-29p is the discrete complex 
Ginzburg-Landau equation. 

In the rest of this subsection, we clarify a part of the phase space for the derived 
system (j2-29p . Using the renormalization equations, we predict the behavior of 
motion in the original system and confirm it numerically. Here we restrict ourselves 
to the conditions Xj+Ar(t) = Xj{t) with N being the number of the oscillators, this 
conditions lead to Ajj^]\j{t) = Aj{t). There is the trivial solution Aj{t) = in Eq. 
(j2-29p . We show that the uniform solution, expressed as Aj{t) = A{t), (for any j), 
can be viewed as one of the local stable manifolds of the fixed point Aj = when 
a certain condition is satisfied. To do this, we study the linear stability for Aj = 
0. Substituting Aj{t) = aj{t), (\aj{t)\ <^ 1) into Eq. ()2-29p we have the linearized 
equation of motion in Fourier space 

dbk{t) -iev r (2'Kk\ 



dt Q 



|_l + e-^'-cos(^^)|6fc(t), (2-30) 



where 



N-l N-l 

^^(^) = ^ E e-^^-^^/%(t), ~a,{t) - ^ E ^''^''"^b,{tl 
^ j=o ^ k=0 

with (A; = 0,...,A^ — 1). From Eq. (j2-30p . we can predict which modes increase or 
decrease in time. The absolute value of bk{t) decreases to zero as time evolves when 

— sin(I2r)cosf — J >0, (2-31) 

and this condition with A; = gives that the uniform solution can be viewed as the 
local stable manifold of Aj = 0. Fig. [2] shows that our analysis for the given system, 
via the renormalization equation, is valid. 

2.5. Spatially extended system (1) 

To show that our renormalization method is also useful in the case that a spa- 
tially extended system modeled by including unperturbed terms with delayed argu- 
ments, we consider the following system, 

du(t,x) TT , . ( , ^-x d'^u(t,x)\ 

^' ^ +—u{t-r,x) =£{au{t,xf + iy „ ^ M , (2-32) 



dt 2r ' \ ' dx"^ 

where a, z^(g M) are parameters, e(€ M) is the small parameter, and r(G M) represents 
the time-delay. Although the system is described by a delay partial differential 
equation, our prescription is not changed. 

First, the naive perturbation solution, x) = u^^> {t,x) +eu^'^\t,x) + ©(e^), 
is obtained by solving the following equations, 

= 0, L,n(i) (t, x) = au(°)=^ (t, x) + t^ ^'"!"^ 1^' , 

du(t,x) -Kuit — r^x) 
L.uit,x) :=4^ + ^^- 
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Fig. 2. The time sequence of the system described by Eas. (l2-26p and (|2-27|) . The values of the 
parameters are J7 = 0.5, a = 1,1^ = 1.01, r = 1 and e — 0.01. The number of the oscillators A'' 
is three, and Xj+N{t) = Xj{t). The initial conditions are Xj{t) = 0.1, pj{t) = 0, {j — 0, 1,2) for 
— r < t < 0, which correspond to the uniform solution. The amplitudes of Xj(t),{j = 0,1,2) 
decrease to zero as time evolves, which we can predict using the condition (|2-3ip with fc = [ 
See text ]. The numerical simulation method is given in the caption to Fig[l] 



The solutions are given by 

1 ~t~ ^ TT 



where A{x){^ C) is an arbitrary differentiable function of x. The naive perturbation 
solution includes the secular term. 

Next, to remove the secular behavior, the renormalized variable A[t,x){^ C) is 
defined as 

x) := A{x) + e-^ f 3a| + i^^^) . 

1 + Z-2 V OX'^ J 

Finally, the renormalization equation, which A{t, x) should satisfy, is derived as 

dAit,x) £ ( ~, ,,2~, , d'^A(t,x) 
^ - ■[3,a\A{t,x)\^A{t,x) + v- ^ ' ' 



dt 1 + if V dx"^ 
This is the well-known complex Ginzburg-Landau equation. 

2.6. Spatially extended system (2) 

To show that the renormalization method is useful in the case that a spatially 
extended system modeled by including perturbation terms with delayed arguments, 
we consider the following system, 

52 A . 3, 



-Q^- ^ + lju{t,x) =£au'^{t-r,x), (2-33) 

where a(G M) is a parameter, r(G M) is the time-delay, and e(G M, |e| <C 1) is the 
small parameter. Along with the procedure for partial differential equations proposed 
in Ref.l6), we can derive the reduced equation. 
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First, the naive perturbation solution u{t,x) = u^'^\t,x) + £u^^\t,x) + ©(e^) is 
obtained by solving the following equations, 

= 0, Lu^^\t,x) = au^°^^{t-r,x), 
Lu(t,x) := {dl - dl^ + l)u{t,x). 

We set n(°) (t, x) as 

n(°)(t,x) = Ae*^'^^-"*) + C.C., 

where A(e C) is the integration constant, and cj^ = A;^ + 1, with A;(g M) being a 
parameter. A secular solution of u^^\t,x) is found to be 

= 3a|yl|^e*'^^XPioi + Poi2;)e^('=^-^*) + C.C., 

where j>io(£ C) and poi(£ C) are parameters. The values of these parameters are 
restricted by the condition 

2{ujpiQ + kpQi) = i. (2-34) 
Next, the renormalized variables A{t,x){^ C) is defined as 

I(t,x) := A + e'ia\A\'^Ae''^''{piQt+poix). (2-35) 
From this definition (j2-35p . we have 

dtA{t, x) = 3ae|I(t, a;)e*'^>io, (2-36) 

d^A{t,x) = 3ae|I(t,x)pI(t,a;)e*'^>oi- (2-37) 

Finally, eliminating pi^ and poi from Eqs. (|2-34p . (j2-36p and (|2-37p . we have the 
following renormalization equation, 

§3. Extended Renormalization Method 

In this section, we propose an extended renormalization method which can lead 
to a reduced equation with delay from a given system with large- or order-one de- 
lay. In this paper, by large-delay we mean that the delayed arguments are of order 
l/e",(Q! > 0) with £ being the small parameter associated with the given weakly 
nonlinear system. We show that our reduced equations are consistent with those 
obtained by the multiple-scale method in Refs.23), 24). 

3.1. Linear system 

The model which we study here is 

"^'^^^^ +a;2x(t) + exft-— ^ =0, (3-1) 



dt^ 

where r/e" represents large-delay with e(G R) being a small parameter, and a(> 0) 
is a parameter. When a = 0, this equation is the same as Eq. ()2Tp . When q = 1 and 
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cv = 1, the system (13-ip was analyzed using the multiple-scale method in Ref23). We 
analyze this system (j3-ip using the extended renormalization method and compare 
the result with that in the previous work. 

First, we obtain the perturbation solution, x{t) = x^^\t) +ex^'^\t) +e'^x^^\t) + 
0{e'^), by solving the following equations 

Lx(0)(t) = 0, (3-2) 
Lx«(t) = -x(°)ft-^Y (3-3) 



it) = -x(i) (t - ^) , Lx{t) := + u.^) x{t). 

In deriving these equations, the magnitude of r/e" in the delayed argument is treated 
as a large value, this treatment corresponds to the nonstandard expansion in the 
previous study. 

The solution of Eq. (|3-2p is 

(t) = ^(0) exp{iujt) + c.c, (3-4) 

where A{0){£ C) represents the contribution to the solution x^'^\t) except for the 
fast motion exp(iu;t). We assume that the solution x^^^ at t — r/e" is expressed as 

=^(^- ^^e-^'^"/^%xpM) +C.C.. 

Here A{—r/e°') represents the contribution to x^'^\t—r/e°'), except for e"*'^^/'^" exp(zu;t). 
This implies that the argument of A{t) is only affected by a large time shift. At this 
stage, we do not know the functional form of A{t). The existence of A{t) is the most 
fundamental assumption in this extended method. It is noted that the equation 
which A{t) should perturbatively satisfy is our extended renormalization equation. 
This extended reduced equation is constructed by removing the secular behavior 
coming from the resonance between the frequency in the operator L and oc exp(iu;t) 
in the forcing terms. When the delay r becomes zero, the extended method corre- 
sponds to the conventional method. Substituting this solution x^^\t) into Eq. (j3-3p . 
we obtain 

Lx«(t) = -^(^- ^^e-^'^^/^%xpM) +C.C.. 
The solution is given by 

xW{t) = —a( - —]e-''^'^/'" exp{iojt) + c.c. 



2uo \ e 
At the next order in e, we obtain 



We observe the secular behavior as we have already seen in the case that the delay 
is not large. 
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Next, we remove the secular behavior. To do this, we define the extended renor- 
mahzed variable. 

This definition is a form of a near-identity transformation at the function ^(0), 
instead of that at the constant A in the conventional renormalization method. From 
this definition of A{t)^ we obtain 

A{t + a) - A{(j) _ ei_^^ _ ^_^^/^. 



a 2uj \ e 



where o"(s M) is a parameter whose value is of smaller than r/e". The inverse of 
Eq. ()3-5p is derived as 

^(0) = A{t) - e^A(t - ^)e--/^'^ + 0{e'). 

Using the above expression, we obtain the following relation 

r 

t — m — 

'-^^^^p^A{t - m^) e— /.^ + 0{e\ (3-7) 

with m(E N). We substitute Eq.(l32]) into Eq. ifM]) an^ take the limit cr ^ 0, 
we obtain our extended renormalization method which A{t) should perturbatively 
satisfy, 

dA{t) i -i^^rle°' 7(4. ^ \ , ,2 ~^ -2iujr/£°' jl 







— m — 


= Aiy 











dt 



Ae-W^"lf t--\+ e2^e-2We" i - 2- V (3-8) 



In Eq. (j3-8p there are delay terms, and this renormalization equation in the case of a = 
1 and w = 1 is equivalent to reduced equations derived in Ref.23), where numerical 
simulation and some analysis have shown that the reduced system reproduces the 
behavior of slow motion in the original system. In the case that a = and r is given 
by Eq. (|2-3|) . we can show that one of the solutions to Eq. (|3-8p up to 0{e'^) is given 
by Eq. (1240]) . 

Here we compare this extended renormalization method with the conventional 
method discussed in ^for this svstem (|3Tp . When we use the conventional method 
we cannot obtain a reduced equation. To see this, we use the conventional method. 
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The naive perturbation solutions are 



3,(i)(t) = iiAe-^-/."e^'^* + c.c. 

+C.C. 



The renormalized variable is defined as 

f2 



^ ^ 2a; \8u;3 \^8a;3 4e" J J 

The renormalization equation up to 0{e) becomes 



dA{t) _ ^i-iujv/s 



dt 2UJ ^ " ^ ^ 

and that up to 0{e'^) becomes 

^ = e±e-^^^/'Ait) + e'( l-V-2W^'^I(t). (3.10) 

Since the magnitudes of the terms calculated as higher-order correction in Eq. (|3T0p 
are 0{e^) and 0(6^^°), this approximation is in contradiction with Eq. (|3-9p except 
for the case of a = 0. When a = 0, Eq. (j3-10p becomes Eq. (|2-8p . and there is no 
contradiction only in the case a = 0. We conclude that, for systems with large-delay, 
the extended renormalization method should be used. 

3.2. Nonlinear system 

We consider a weakly nonlinear system with large-delay which appears in optics. 
In Ref.24), they have analyzed the system with optoelectronic feedback, and the 
system is described as 

-y{s)-eMs)(i + T^y(^))^+^'cli + y(s ® 



dy{s 



ds 



{l + y{s))x{s), 



where s is the scaled time C, P, 0{€ M) are parameters, and e{£ M) is the small 
parameter. The solution which we focus on is the small amplitude regime, described 
by the following assumption 

x{s) = ex(i)(s) + e^x^^\s) + e^x^^\s) + O(e^), 
y{s) = ey^'\s) + e^y^^\s) + e^y^^\s) + 0{e^). 

We construct the reduced equation using our extended method, and compare the 
result with that reported in Ref.24). 
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First the naive perturbation problems are 



-y^-'>is)-x^'>is) + Cy^'>is-9), 



ds 

where = Qje^ . The solutions x^^''(s),x*-^^(s) and x^'^^(s) are given by 
x(i)(s) = ^(0)e^^ + c.c, = Z!^(0)2e2^^ + c.c, 

x(3)(.) = I AC7^(0) - ^|^(0)|M(0) -^(0) 



—iCA{—Q)j^'^ + c.c. + higher harmonics. 
The definition of the renormalized variable is 

A{t) := ^(0) + e^f i^CA{^) - ^|^(0)|M(0) - ^(0) - iCA{-Q)^ . (3-11) 

The renormalization equation is derived from Eq. (j3-lip as 

= -( iCA{s) - -\A{s)fA{s) - A{s) - iCA{s - 6)] . (3-12) 
as 2 \ 3 / 

The renormalization equation (j3-12p is the reduced equation derived in Ref.24). Some 
analytical analyzes in Ref.24) have shown where bifurcation points are. Again, we 
confirm that our extended method gives the same results given by the use of the 
multiple-scale method. 

§4. Discussion and Conclusions 

In this section, we discuss the features of both our conventional and extended 
renormalization methods, and then we conclude our study. 

The prescription of the conventional renormalization method for systems with 
order-one delay is not different from ones without delay. The standard prescrip- 
tion leads to the reduced equation from a given weakly nonlinear system. The 
conventional method removes the secular terms from naive perturbation series by 
accounting for their effect with renormalized variables. Derived reduced systems are 
always ones without delay. Being without delay in a reduced equation means that 
the dimension of phase space for the original system can be reduced perturbatively. 
This reduction provides us the approximate structure of phase space. Systems to 
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which we can apply this method are weakly nonlinear ones with order-one delay. In 
this sense, the conventional method is restricted. 

The prescription of our proposed extended rcnormalization method also removes 
the sccularity. The basic assumption for this method is that we can introduce an 
unknown function contributing the naive perturbation solution, instead of the inte- 
gration constant in the use of the conventional method. Although a rigorous math- 
ematical meaning of the extended method has not yet given in this paper, we have 
checked the validity of our method through various examples. Derived reduced sys- 
tems using this extended method are always ones with delay. This means that the 
dimensions of phase space for both the original and reduce systems are high. The 
advantage of our reduction method is that a steady state in the reduced system cor- 
responds to a periodic one in the given system, which provides us some bifurcation 
analyzes. Using extended method, we can deal with systems whose delay time is of 
order 1/e" , (a > 0) where e is the small parameter appearing in the original system 
under study. 

For both the rcnormalization methods, terms in the reduced equation arise from 
secular terms appearing in the naive perturbation analysis. This implies that higher 
harmonics in the naive perturbation analysis does not contribute to the reduced 
equation in the first order approximation. In this sense, the reduced equation can be 
obtained from a wide class including the original system. Compared to the multiple- 
scale method, our methods do not need scaled variables. While in the course of the 
derivation of a rediiccd system Tising oTir methods, we need the analytical expressions 
of the naive pcrtTirbation solutions so that we define the rcnormalized variables. Al- 
though the procedures of our methods are systematic, the application of our methods 
are restricted by this disadvantage. 

In this paper, we have shown that the rcnormalization method can be extended 
to a tool to study systems with delay, and that the method gives reduced systems 
successfully. Combining the previous studies of the rcnormalization method with the 
present study, we expect that our rcnormalization method includes all the asymp- 
totic analyzes. Furthermore, we believe that the application of the renormalization 
method can help elucidate the behavior of time-delayed systems in a non-chaotic 
regime. 

Acknowledgements 

The author would like to thank the members of NTT Communication Science 
Laboratories for their continual encouragement. 

References 

1) R. Lang and K. Kobayashi, IEEE J. Quantum Electron. 16 (1980), 347. 

2) D. Pieroux and P. Mandel, Phys. Rev. E 67 (2003), 056213. 

3) A. Uchida, N. Shibasaki, S. Nogawa and S. Yoshimori, Phys. Rev. E 69 (2004), 056201. 

4) G.F. NcwoU, Opr. Res. 9, (1961) 209; G. Orosz, R.E. Wilson and B.Krauskopf, Phys. Rev. 
E, (2004) 026207. 

5) M. C. Mackey and L. Glass, science 197 (1977), 287. 

6) M. C. Mackey, J. Econ. Theory 48 (1989), 497. 



Renormalization Reductions for Systems with Delay 



17 



7) K. Gopalsamy, Stability and Oscillations in Delay Differential Equations of Population 
Dynamics (Kluwer Academic, 1992). 

8) M. A. Johnson and F. C. Moon, Int. J. Bifurcation and Chaos 9 (1999), 49. 

9) L. Y. Chen, N. Goldenfeld and Y. Oono, Phys. Rev. Lett. 73 (1994), 1311; Phys. Rev. E 
54 (1996), 376. 

10) Y. Oono, Int. J. Mod. Phys. B 14 (2000), 1327. 

11) T. Kunihiro, Prog. Theor. Phys. 94 (1995), 503; ibid 97, (1997) 179, T. Kunihiro and J. 
Matsukidaira, Phys. Rev. E 57 (1998), 4817. 

12) S-I. Ei, K. Fujii and T. Kunihiro, Annals. Phys. 280 (2000), 236. 

13) B. Mudavanhu and R.E. O'Malley, Jr., Stud. App. Math. 107 (2001), 63. 

14) K. Nozaki and Y. Oono, Phys. Rev. E 63 (2001), 046101. 

15) M. Iwasa and K. Nozaki, Prog. Theor. Phys. 116 (2006), 605. 

16) S. Goto, Y. Masutomi and K. Nozaki, Prog. Theor. Phys. 102 (1999), 471. 

17) S. Tzcnov, Comtemporay accelerator physics (World Scientific, 2004). 

18) Y. Nambu, Phys. Rev. D 65 (2002), 104013; Y. Nambu and Y.Y. Yamaguchi, Phys. Rev. 
D 60 (1999), 104011. 

19) M. Frasca, Phys. Rev. A 56 (1997), 1548. 

20) Y. Masutomi, and K. Nozaki, Physica D 151 (2001), 44; S. Kawaguchi, Prog. Theor. Phys. 
113 (2005), 687; I. Rozhkov and E. Barkai, J. Chem. Phys. 123 (2005), 074703. 

21) S. Goto and K. Nozaki, Prog. Theor. Phys. 105 (2001), 99; J. Phys. Soc. Jpn. 70 (2001), 
49 ; S. Goto, K. Nozaki and H. Yamada, Prog. Theor. Phys. 107 (2002), 637; S. I. Tzenov 
and R. C. Davidson, New J. Phys., 5 (2003), 67; S. Goto and K. Nozaki, Physica D, 194 
(2004), 175; T. Maruno, S. Goto and K. Nozaki, Prog. Theor. Phys. Ill (2004), 463; S. 
Goto, Prog. Theor. Phys. 115 (2006), 251. 

22) P.M. Alsing, V. Kovanis, A. Gavrielidcs and T. Erncux, Phys. Rev. A 53 (1996), 4429. 

23) S.L. Das and A. Chatterjee, Nonlinear Dynamics 39 (2005), 375. 

24) D. Pieroux, T. Erneux, A. Gavriehdes and V. Kovanis, SIAM J. Appl. Math. 61 (2000), 
966. 



